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Lepton flavor changing in neutrinoless r decays. 
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Abstract- Neutrino oscillations, as recently reported by the Super-Kamiokande collaboration, 

imply that lepton numbers could be violated, and — > ^ + £ + + l~ , — » ^ + p° are some 

Q H ' typical examples. We point out that in these neutrinoless modes, the GIM cancelation is much milder 
i \ 

^h. with only a logarithmic behavior log(m :7 /mfc) where rrij t k are the neutrino masses. This is in sharp 
contrast with the vanishingly small amplitude — > +7 strongly suppressed by the quadratic power 
(rrij — mfy/M^j. In comparison with the hopelessly small branching ratio B(r =t — ► + 7) ~ 10 -40 , 
the B(r =t — > ^ +£ + +£~) could be larger than 10 -14 . The latter mode, if measurable, could give one 

_ 1 

more constraint to the lepton mixing angle sin 26 jk and the neutrino mass ratio rrij/mk, and therefore 
is complementary to neutrino oscillation experiments. 
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Evidence for the transmutation between the two neutrino species <-> v T is recently reported by 
the Super-Kamiokande collaboration^ . As a consequence, neutrinos could have nondegenerate tiny 
masses and lepton numbers would no longer be conserved. Hence, besides the well known neutrino 
oscillation phenomena, neutrinoless r decays, such as — > ^ + 7, — ► fj^ + £ + + £~ and — > 
^ + p° could occur. 

The interest of the r^ 1 — > ^ + £ + +£" and/or r 1 * 1 — > + p° modes is twofold. First, contrarily to 
the radiative case — ► /r t + 7 which is damped by a vanishingly quadratic power 2 , the typical lepton 
flavor changing amplitude — > + £ + + t~ (or — > ^ + p°) is only suppressed by a smooth 
logarithmic term. Second, these r decay modes, if measurable, are complementary to the neutrino 
oscillation experiments. They could give - besides the lepton mixing angle 9ij - the ratio m 2 /m|, 
whereas neutrino oscillations give the difference \m 2 — m\\. By combining them, the absolute value of 
rrij could be determined in principle. 

Similarly to the CKM flavor mixing in the quark sector, let us assume that the neutrino gauge- 
interaction eigenstates v e , and v T are linear combinations of the three neutrino mass eigenstates 
v\, v 2 and z> 3 of nonzero and nondegenerate masses m-i, m 2 and m 3 respectively. Thus 



(1) 



where the 3x3 matrix U\ cp is unitarity. Neutrino oscillation measurements give constraints usually 
plotted in the (sin 29, Am 2 = \m 2 — m 2 \) plane, where 9 is one of the three Euler angles in the rotation 
matrix U\ ep . 

The weak interaction effective Lagrangian for charged current of leptons can be written as 
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where the charged current L\ is 



j'=i 



Here i stands for e~ , yT , t~ and uj (with j = 1,2,3) are the three neutrino mass eigenstates. For any 
fixed £, one has J2j \ Uej\ 2 = 1- For instance the v^, operationally defined to be the invisible particle 
missing in the ir + — > /U + + v^, is initially a superposition of v\, V2 and ^3, in the same way as the K° 
meson produced by strong interaction, say by ir~ + p —> K° + A, is initially a superposition of the 
mass eigenstates Kl and Kg with masses mi, / ms- The nondegenerate masses are the origin of the 
oscillation phenomena of both neutrinos and neutral K mesons. 
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In the most general renormalizable gauge, at one loop level to order g A - where g = e/ sin#w 
is the weak interaction coupling constant - there are in all eighteen Feynman diagrams contributing 
to the neutrinoless decays — ► /z^ + i + + l~ or — ► /x^ + p°; these modes are mediated by the 
Z and the photon, ten diagrams for the virtual Z, and eight for the virtual 7. Three of them are 
depicted in Figs. 1-3. The fifteen others, not shown here, are similar to Figs. 1-3 in which the internal 
W 1 * 1 in loops are replaced in all possible ways by the "would be" Goldstone bosons $ ± , those absorbed 
by the gauge bosons to render them massive by the Higgs mechanism. The intermediate virtual 
photon is absent in Fig.l and in Fig. Ibis which is the analogue (not depicted here) of Fig.l with 
replaced by < J> ± . The contributions of the mediated neutral Higgs boson H° are negligibly small for 
both reasons (its mass and its couplings with the lepton £ pair or the up down quarks of the p°) and 
can be discarded. 

A careful examination of these eighteen diagrams shows that only Fig.l and Fig. Ibis provide 
the logarithmic behavior log(mj/M^), while the contributions of all other 16 diagrams are power 
suppressed as (mj/M^), {m 2 /M^j) x log(mj/M^) and therefore vanishingly small. The principal 
reason for the appearance of the logarithmic log (nij/M^) term is that we are dealing in Fig.l and 
Fig. Ibis with two propagators of nearly massless fermions for which if the momentum transfer q 2 is 
much smaller than and consequently neglected, infrared divergences appear when mass of the 
internal fermion goes to zero 3 ' 4 . We emphasize that Fig.l and Fig. Ibis are the only ones that contain 
an infrared divergence log(m|/M^), this fact has been noticed longtime ago in different contexts, for 
instance in the computation of the slope of the neutrino electromagnetic form factor 3 , and the s-d-7 
induced coupling 4 . Note however that compared to Fig.l, the contribution of Fig. Ibis is damped by 
an additional Mm/M^ factor because of the $-fermion couplings, where M and m are respectively 
the r lepton and muon masses. So actually only Fig.l dominates. 

Due to the unitarity of the U\ ep reflecting the GIM cancelation mechanism, the divergence as well 
as the mj-independent finite part of the loop integral do not contribute to the decay amplitude because 
they are multiplied by ^2j(U*jU T j) = when we sum over all the three neutrino contributions. Only 
the mj-dependent finite part of the loop integral is relevant. This point is crucial, implying that we 
cannot neglect rrij no matter how small rrij is, otherwise we would get identically zero result after the 
summation over the neutrino species j. 

Our first task is to show that Fig.l and Fig. Ibis actually give rises to the logarithmic log(mj/M^) 
term. This term could be equally guessed by approximating the W propagator with i/M^, the W 
mass plays the role of the loop integral momentum cutoff. Hence Fig. 1 looks like the familiar fermionic 
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loop of the gauge boson self energy, or vacuum polarization. When q 2 = (q being the four-momentum 
of the external gauge boson), the standard log(m 2 /^ 2 ) appears 5 . The following calculation of the 
diagram of Fig.l confirms this expectation. 

Let us write the one-loop effective T-//-Z transition of Fig.l as u(p)Tj(q 2 )u(P)e\(q), where P, p 
and q = P — p are respectively the four-momentum of the r lepton, muon and virtual boson Z. Thus 

2 ' d 4 fc 
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The £ dependence in the W propagator of Fig.l is canceled by the £ dependence of Fig. Ibis (where $ 
replaces W), so for simplicity, we take the £ = 1 Feynman-'t Hooft gauge at the outset. 

Inserting Tj(q 2 ) inside u(p) and u(P), making use of Dirac equations for these spinors and adopting 
the standard Feynman paramerization for the denominator in Y*(k, q), we get after the k integration 
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The form factors a, ■ ■ ■ , f are 



a = Vj log 



2l 

M 2 



+ M 2 x(x + y) +m 2 x(l - y) - q 2 (x + y)(l - y) , b = Mmx , c = —M 2 x(x + y) 



d = Mmx(l + y) , e = —M 2 (x + y)(x + 2y - 2) , / = Mm [x(l + y) - 2y(l - y)] . 



(6) 



The log(Pj/M^) in (6) is the finite part extracted from V~ £ T(e) « [1 - elog(P i /M 2 v )]r(e) where 
e = 2 — (n/2) is the n-dimensional regularization parameter used to handle the ultraviolet divergence 
of (2). As discussed, the j-dependence in T>j(q 2 ) as given by (5) is crucial to get nonzero result. On 
the other hand, since the j-independent terms q 2 y(l — x — y) < q 2 /4 < (M — m) 2 /4 <C is much 
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smaller than M^, it is useful to make an expansion of T x (q 2 ) in power of rj = q 2 /M^ by writing 
Vj(q 2 ) = Vj(0) +0(rj), where 

Vj(0) = x [M$f -m 2 - M 2 {1 -x-y)- m 2 y] + m 2 . (7) 

This rj expansion simplifies the x, y integrations and shows us that the relevent j-dependent part of 
T x (q 2 ) has the following general form (the Lorentz index A is omitted for simplicity) : 

Tj(q 2 ) = (Ar, + B)F(r,,m j ) + G(m j ), (8) 

where A, B are constant, and most importantly we note that F(0,mj) = log^-, with Sj = m 2 /M^. 
The second term G{mj) only contains Sj, Sj log5j and S 2 log5j terms which are negligible for rnj — ► 0. 

The general structure in (8) emerges after the x,y integrations of (4), using (7). This dominant 
log Sj manifests itself from the lower limit x = of the / / dxdy [1/Pj(0)] integration. Note that 
/ J dxdy [1/Pj(0)] gives log<5j, while / / dxdy [x n /Vj(0)} gives SJ log<5j. Therefore, if there exists 
terms independent of x in the numerator M x (q 2 ), then these x-independent terms will lead to log Sj. 
There are actually three x-independent terms in Af x (q 2 ) - those proportional to y(l — y) - localized 
in the form factors a, e and / which give rise the log<5j after the / J dxdy [M x (q 2 )/Vj(0)] integration. 
The constant A in (8) comes from the form factor a and the constant B from the form factors e, /. 



The x-dependent terms in N (q ) yield Sj, Sj log Sj and S 2 log Sj terms, they are grouped into G( 



The presence of this infrared divergence log Sj from diagrams similar to Fig.l has been noticed in 
the literature 3 ' 4 . We remark that the q 2 y(l — y) term (and not the log(T>j(q 2 )/M^)) in the form factor 
a of (5) which yields the important log Sj. The form factors e and / on the other hand, because of their 
q x operator, contribute negligibly when contracted to the lepton pair current. The dominant 

term of T x {q 2 ) is found to be 



6M 2 , 
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From (9), the r — ► \i + i + + l~ decay amplitude is 
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C/ TJ .C/*.log-^-] 7 A (l-75). (9) 




L x = u(p) 7 A (l - j 5 )u(P) , i\ = u(k_)j x (gv ~ 759A)v(k + ) , g v = — + 2 sin 2 8 W , g A = 

For A(t — > (i + p°), we simply replace i\ by m p f p e\ where f p ps 150 MeV is the decay constant of the 

p , extracted from p° — > e + + e~, and £>, is the p polarization vector. It remains to evaluate 

3 / m 2 \ 3 / 2\ 

B = E ^ log ( ^ J = E UrkU* k log ( -J ) (12) 



where we have used J2j U^jU T j = to get rid of the first V\ mixing parameter C/'* 1 C/ T i. The factor B 
in (12) which represents the soft GIM cancelation tells us that when rrij are degenerate, the lepton 
flavor mixing does not occur and the neutrinoless-r decays identically vanish. 

To estimate B, let us assume^ the following form of the Ui ep , neglecting possible CP violation in 
the lepton sector: 
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The mixing angle 023 ~ 45° is suggested by the Super-Kamiokande data and the 0i3 « 0° comes 
from the Chooz data^ 6 ) which give 0i3 < 13°, whereas 9u being arbitrary. Thus, 

B = cos 2 12 log — + sin 2 12 log — . (14) 
m 3 m 3 

Although 0i2 is likely small « 0°, however the maximal mixing 0i 2 ps 45° may be also possible allowing 
v e v a (as suggested by the LSND experiment). Taking 0i2 in the range 0°-45°, and using^ 1 ' 6 ) 
Am^ = \ m 1 ~~ m il = 2 x 10~ 3 eV 2 , m 3 w 5 x 10~ 2 eV, then \B\ 2 is of order of unity [\B\ 2 ~ 0(1)}, it 
could be bigger if m\ or 777-2 are exponentially smaller than 7773. We get the branching ratio 

B(r -»• fi + £ + + r) > 10~ 14 . (15) 

Although (15) is still very small, however compared to the radiative case r — > \i + 7, it represents 
a spectacular enhancement through the mild GIM cancelation. Finally, we remark that for large 
q 2 > M(y, for instance in /x^ + r T , we have the same ten diagrams with a real Z boson. However, 
obviously the q 2 = M| expansion does not make sense, and F(M| , mj) in (8) cannot be approximated 
by F(0, rrij) = \og5j. Only for small q 2 <C that the use of the dominant term F(0, rrij) can be 
justified. Therefore the low energy e + + e~ — > /j 1 * 1 + reaction might be also worth to investigate. 
However compared to the one-photon exchange e + + e~ — > /^ + + cross-section, the e + + e~ — ► 
/X"'" + r T cross-section is damped, besides the coefficient £>a em /247r sin 2 W squared, by an additional 
s 2 1 (Mf — s) 2 multiplicative factor due to the Z propagator. 

Independently of the precise numerical value of B, the neutrinoless decay modes r — > /j + p° 
and/or r — > /J + £ + + £~ are interesting on their own right for two reasons. First, although being 
higher order loop effect, the branching ratios are not desperately small due to the smoothly logarithmic 
GIM suppression. Second, while neutrino oscillations only provide the mass difference Am^, lepton 
flavor changing r decays give the ratio m 2 jm\. When combining these two processes, the absolute 
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value of the neutrino mass rrij may be obtained. Both reactions are mutually complementary in the 
determination of the neutrino masses and the lepton mixing angles. 

Lepton flavor changing processes in the seesaw-type neutrino models are also discussed in 7 . 
Acknowledgments : Helpful discussions with John Iliopoulos are gratefully acknowledged. 
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Figure Caption : 

Figures 1-3 : One-loop lepton flavor changing r — > fj, + £ + + l~ . 
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